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Abstract  -  Ground  targets  are  constrained  on  the  Earth 
with  their  velocity  vector  direction  aligned  mostly  along 
the  body  longitudinal  axis.  The  pose  angle  therefore 
carries  kinematic  information  useful  for  tracking 
maneuvering  targets.  For  target  identification  (ID), 
range  profiles  obtained  by  a  high  range  resolution  (HRR) 
radar  are  compared  with  reference  templates  in  pose 
angle  per  target  class,  thus  producing  pose  angle 
estimates.  In  this  paper,  we  present  a  method  for 
measuring  the  pose  angle  of  a  maneuvering  target  by 
first  converting  the  matching  scores  of  a  classifier  into  a 
“likelihood”  function  and  then  updating  the  likelihood 
function  probabilistically.  By  accumulating  the 
likelihoods  of  possible  poses  over  time,  it  enables  the  pose 
angular  tracking  as  the  target  is  undertaking  turn 
maneuvers.  Simulation  results  for  pose  angular  tracking 
are  presented  wherein  range  profiles  are  generated  from 
RF  signatures  of  moving  targets. 
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1  Introduction 

Compared  to  conventional  tracking  with  post-detection 
observables  such  as  range,  range  rate,  and  bearings, 
feature-aided  tracking  (FAT)  works  on  low-level 
measurements.  When  compared  with  reference  templates 
in  a  database,  a  successful  matching  provides  target  type 
and  viewing  angles  (or  pose  angles)  among  others.  It  thus 
offers  an  extended  target  state  including  not  only  the 
kinematic  variables  such  as  position,  velocity,  and  possibly 
acceleration  but  also  the  target’s  orientation  relative  to  its 
environment  (terrain  and  road). 

For  an  HRR  radar  [14],  range  profile  is  a  one-dimensional 
(ID)  measurement  of  target  radar  reflectivity  along  the 
radar  to  target  line  of  sight  (LOS)  vector  where  the 
amplitudes  are  statistical  features.  This  look  vector,  when 
expressed  in  the  target  body  frame  in  terms  of  the  aspect 
and  depression  angles,  is  called  a  “pose.”  For  practical 
reasons,  a  target  is  pre-sampled  into  a  template  library  of 
range  profiles  at  discrete  poses.  A  successful  template 
matching  therefore  identifies  the  target  type  and  at  the 
same  time  produces  the  pose  at  which  the  range  profile  is 
generated. 

HRR  range  profiles  have  long  been  used  for  target 
identification  (ID)  or  fingerprinting  [8,  9,  13,  15].  It  has 
also  been  used  in  data  association  to  improve  track 
continuity  [3,  6].  Recently,  the  pose  angular  measurements 
as  a  by-product  of  target  ID  have  been  used  to  assist  target 


Erik  Blasch 

Air  Force  Research  Lab 
WPAFB,  OH  45433 
erik.  blasch  @  wpafb .  afb .  mil 

tracking  [16].  It  exploits  the  fact  that  ground  targets  are 
constrained  to  move  on  the  Earth  surface  and  their  velocity 
vector  direction  is  aligned  mostly  along  the  body 
longitudinal  axis.  The  pose  angle  therefore  carries 
kinematic  information  which  is  particularly  helpful  and 
fast  about  target  maneuvers. 

The  angular  pose  of  a  target  is  typically  estimated  in  two 
ways.  One  way  is  to  consider  the  target  pose  as  an 
intermediate  result  or  a  by-product  (a  nuisance  parameter) 
of  the  target  ID  process  as  described  above.  A  template 
must  have  about  the  same  pose  angles  (azimuth  and 
elevation,  aspect  and  depression)  as  the  target  in  order  to 
reach  a  successful  classification.  Target  pose  is  used  as  a 
search  parameter  in  the  target  ID  process. 

In  the  second  way,  target  pose  is  estimated  explicitly  for  a 
target  with  known  ID.  For  example,  neural  networks  were 
used  for  pose  estimation  with  for  HRR  and  SAR  images 
[7,  11]  where  the  weight  vector  is  obtained  through 
training. 

It  is  well  known  that  ID  HRR  signatures  are  subject  to 
high  variability  due  to  scintillation  effects  (speckles).  That 
is,  multiple  scatterers  fallen  in  a  single  range  bin  interact 
either  constructively  or  destructively  when  the  aspect  angle 
is  changed  slightly.  This  is  the  basis  for  rather  sensitive 
pose  estimation.  In  conventional  classifier  designs,  the 
sensitivity  to  pose  angles  is  minimized  so  that  few 
templates  are  required  to  discriminate  between  different 
classes  of  targets.  However,  in  classifier-based  pose 
angular  estimation,  we  would  like  to  maximize  the 
sensitivity  of  a  classifier  to  pose  angles. 

In  this  paper,  a  particular  classifier,  namely,  the  MTE 
algorithm  [4],  is  used  for  pose  estimation.  Due  to  the 
discrete  nature  of  reference  templates  and  errors  in  range 
profiles,  the  pose  angular  estimate  is  not  single-valued  but 
rather  is  a  distribution  of  matching  scores  at  the  output  of  a 
classifier.  In  this  paper,  we  present  a  method  for  measuring 
the  pose  angle  of  a  maneuvering  target  by  first  converting 
the  matching  scores  of  a  classifier  into  a  “likelihood” 
function  and  then  updating  the  likelihood  function 
probabilistically.  By  accumulating  the  likelihoods  of 
possible  poses  over  time,  it  enables  the  pose  angular 
tracking  as  the  target  turns. 

The  paper  is  organized  as  follows.  In  Section  2,  the  MTE 
classification  algorithm  is  first  described.  In  Section  3, 
pose  angular  estimation  with  the  MTE  classifier  is  detailed 
together  with  simulation  results  presented  to  illustrate  the 
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concept  and  performance.  Section  4  concludes  the  paper 
with  future  work  outlined. 

2  MTE  Method  for  Range  Profile  Comparison 

There  are  different  classifiers  available  for  matching  a 
measured  HRR  range  profile  with  reference  range  profiles 
in  a  template  database.  One  method  is  the  linear  regression 
matching  algorithm  [13].  Another  method  that  is  used  in 
this  paper  is  adopted  from  the  DARPA  MTE  program  [4]. 
The  MTE  algorithm  is  a  minimum  squared  error  (MSE)- 
based  algorithm.  It  has  the  following  features.  It  exhibits  a 
matched-filter  behavior,  compensates  for  differences  in 
profile  gain,  searches  over  range  shifts  to  find  the  best 
profile  alignment  and  includes  weighting  to  increase 
importance  of  high-amplitude  peaks. 

The  inputs  to  the  algorithm  include  a  “Target”  HRR  profile 
in  dB,  P,  with  a  profile  mask,  mp ,  and  a  “Template”  HRR 
profile  in  dB,  T ,  with  a  profile  mask,  mt.  It  further  includes 
an  estimated  mean  clutter  value  for  the  target  profile,  c,  a 
number  of  range  shifts  over  which  to  search,  N,  and  a 
weighting  parameter,  w. 

If  the  template  HRR  profile  and  target  HRR  profile  are 
given  in  amplitude,  p  and  t ,  they  will  first  be  converted  to 
the  HRR  profiles  in  dB,  P  and  T,  respectively.  Similarly, 
the  estimated  mean  clutter  value  for  the  target  profile  in 
dB,  C,  is  also  calculated  from  amplitude  c,  which  can  be 
estimated  from  the  target  profile  P  as  the  median  value  of 
the  off-mask  pixels.  The  calculations  are  given  by: 

Template  profile:  p  =  io^/20and  p  =  20 lo g10(p)  (la) 

Target  profile:  t_  =  10-/2°  and  T  =  201og10(f)  (lb) 

Mean  clutter:  c  =  10c/2°  and  c  =  201og10(c)  (lc) 

The  MTE  algorithm  consists  of  the  following  steps.  In  the 
first  step,  the  template  profile  in  T  at  off-mask  range  bins 
is  set  to  the  mean  clutter  C.  The  second  step  is  to  align  the 
profiles  in  range  using  the  profile  masks.  The  two  profiles 
can  be  aligned  either  at  the  leading  edge  of  the  masks  or  at 
their  centroids.  If  the  two  profiles  are  of  different  lengths, 
the  unaligned  parts  are  filled  with  the  mean  clutter,  C.  The 
next  step  is  to  calculate  a  complementary  weight  vector  for 
the  target  profile  by: 

exp(-10(--C)/2°  /  w)  =  exp (-p/cw),  P>C,  on  -  mask 
exp(-l/w),  P<C,  on -mask 
C,  off  -  mask 

From  (2),  the  complementary  weighting  Wp  varies 
according  to  the  profile  amplitude  as  compared  to  the 
mean  clutter.  Off  mask,  the  weight  is  assigned  a  constant 
value  and  on  mask,  when  a  target  profile  value  is  equal  to 
or  smaller  than  the  mean  clutter,  the  weight  is  also 
constant.  When  a  target  profile  value  is  greater  than  the 
mean  clutter,  its  ratio  with  the  mean  clutter  is  taken  as  the 
exponent  of  a  negative  exponential,  and  the  resulting 
complementary  weight  is  rather  small. 

The  next  step  is  to  shift  the  template  profile,  adjust  its  gain, 
and  calculate  its  weight  prior  to  comparing  it  with  the 


target  profile.  There  will  be  2/V+l  shifts  from  -N  to  N.  This 
is  done  by  the  following  loop: 

For  s  =  -N  to  N: 

•  Shift  T  in  range  by  s  range  bins. 

•  Compute  the  average  difference  between  T  and  P  over 
the  intersection  of  their  masks  for  value  above  the  mean 
clutter  C  to  obtain  the  optimal  gain  adjustment. 

•  Adjust  T  by  biasing  it  with  the  average  difference 
(addition  in  dB  is  equivalent  to  scaling  in  amplitude  in 
the  optimal  gain). 

•  Compute  the  weighting  Wt  using  the  same  formula  as  in 
(2)  for  the  gain  adjusted  template  T. 

•  Compute  the  combined  weight  as: 

W=l-W,°Wp  (3) 

where  0  stands  for  element  to  element  multiplication  of 
two  vectors. 


•  Compute  the  MTE  score  for  this  shift  as: 


Score(s)  = 


^Wo(P-T)2 


(4) 


where  the  summations  are  carried  out  over  the  vector 
elements  and  the  summation  in  the  denominator  is  over 
those  bins  with  value  over  c. 


End 

The  final  MTE  score  is  computed  as: 

MTE  =  min  Score(s)  (5) 

Since  the  score  is  marked  by  the  total  distance  between  the 
two  profiles,  more  importance  (large  weight)  is  placed  on 
characteristic  bins.  In  other  words,  the  profile  values  at 
those  bins  have  to  be  close,  otherwise,  they  weight  up  the 
cumulative  differences  between  the  two  profiles. 

3  Pose  Angular  Estimation  with  MTE 

The  MTE  classifier  output  is  also  called  the  MTE 
discriminant  or  discriminant  for  short  and  will  be  analyzed 
in  this  section.  Two  methods  for  angular  estimation  are 
then  described. 


3.1  MTE  Classifier  Discriminant  Analysis 

The  range  profiles  of  two  targets,  i.e.,  a  generic  tank  and  a 
vehicle  transport,  are  generated  with  the  simulation 
environment  shown  in  Fig.  1.  It  consists  of  the  Simulation 
Tool  for  Advanced  Radar  Systems  (STARS),  developed  by 
ATK  Mission  Research  (Dayton,  OH)  [12]  and  the  Radar 
Signature  Predictor  (SigPred),  developed  by  General 
Dynamics  Advanced  Information  Systems  (Ann  Arbor, 
MI)  under  the  Air  Force  program  Feature- Aided  Tracking 
of  Stop-move  Objects  (FATSO)  [10].  More  information 
about  the  simulation  environment  can  be  found  in  [18,  19]. 

In  our  simulation,  the  radar  platform  is  at  an  altitude  of 
5000  meters  at  a  constant  speed  of  100  m/s.  The  radar 
sensor  is  side  looking  at  a  fixed  depression  angle  of  12 
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degrees.  The  center  frequency  is  10  GHz  with  a  bandwidth 
of  600  MHz  and  the  pulse  repetition  frequency  is  2000  Hz. 

The  discriminant  provides  a  measure  of  the  “distance” 
between  the  test  profile  and  the  template  profile.  The 
smaller  the  discriminant  output  or  the  score,  the  closer  the 
two  profiles  look  alike.  Figs.  2  and  3  show  the 
discriminator  outputs  where  the  test  profile  and  the 
reference  templates  are  for  the  same  target  (i.e.,  auto 
correlation)  for  a  generic  tank  and  a  vehicle  transport, 
respectively.  The  cross  correlation  between  the  two  is 
shown  in  Fig.  4.  To  generate  these  plots,  the  target  is 
sampled  in  azimuth  in  a  one-degree  increment  over  360 
degrees.  In  the  auto  correlation  cases,  the  test  azimuth 
angles  are  off  the  template  azimuth  angles  by  0.2  degrees. 

From  Fig.  2,  the  high  correlation  exists  for  azimuth  angles 
6  and  360-  6,  a  longitudinal  symmetric.  This  shows  up  as 
the  diagonal  and  anti-diagonal  components,  a  big  X.  We 
can  also  see  some  minor  correlations  at  angles  180-0  and 
180+0,  which  form  a  big  O.  The  minor  correlation  is  less 
obvious  as  O  in  Fig.  3  but  rather  more  like  small  x.  There 
are  also  weak  correlation  features  around  the  plot. 

There  are  no  distinct  features  in  Fig.  4  for  cross  correlation 
between  generic  tank  and  vehicle  transport  except  wide 
horizontal  strips  at  90  and  270  degrees,  indicating  the  two 
targets  are  difficult  to  discriminate  when  viewed  from 
those  azimuth  angles  (broad  side).  Overall,  the  vehicle 
transport  seems  to  be  easier  than  the  generic  tank  for 
classification  and  pose  estimation. 

The  discriminant  outputs  at  six  consecutive  azimuth  angles 
are  shown  in  Fig.  5  for  the  vehicle  transport  where  the  true 
target  azimuth  angle  varies  from  45  to  50  degrees.  The 
target  range  profile  at  each  azimuth  is  compared  with  the 
set  of  templates  from  0  to  359  degrees  but  only  the  first 
180  results  are  plotted.  The  discriminant  values  (the 
distance  between  the  test  profile  and  the  template  file)  are 
the  lowest  around  the  true  azimuth  angle.  There  is  a 
secondary  dip  around  95  degrees  due  to  minor  symmetry. 

Fig.  6  is  the  blow-up  of  Fig.  5  around  the  dips.  When  the 
index  of  the  lowest  dip  is  used  as  the  pose  estimates  (see 
Section  3.2  for  further  discussion),  the  angular  estimation 
at  some  azimuth  will  develop  errors.  Also  note  that  the 
discrimination  curves  do  not  have  the  same  slopes,  some 
are  steeper  than  others.  Besides,  the  left  and  right  slopes 
are  not  symmetric.  Those  azimuth  angles  that  produce 
steep  discrimination  curves  can  be  estimated  to  a  better 
accuracy. 

In  summary,  when  the  MTE  algorithm  is  used  for  pose 
estimation  via  search,  the  discriminant  output  exhibits 
certain  azimuth  symmetry.  This  includes  left  vs.  right 
(symmetry  about  longitudinal  axis),  front  vs.  back 
(symmetry  about  the  perpendicular  axis),  and  some  local 
symmetries.  In  addition,  the  discriminant  output  is  volatile 
(large  variations  in  the  score  values  as  a  function  of 
azimuth  angle)  and  peaks  do  not  have  regular  shapes. 

The  ambiguity  due  to  symmetry  may  be  solved  by  the 
target  velocity  vector,  road  direction,  and  sensor  to  target 


geometry.  The  search  interval  in  azimuth  may  be  restricted 
to  a  small  sector.  However,  there  is  a  need  to  consider 
target  maneuver,  which  can  change  heading  rather  quickly 
particularly  over  a  long  sampling  interval.  Several  pose 
estimation  algorithms  are  presented  below. 

3.2  Discrete  Pose  Estimate  at  Minimum  Distance 

The  discriminant  typically  does  not  show  a  clear  winner 
over  a  search  interval.  Therefore,  there  are  several  possible 
ways  to  obtain  an  estimate: 

•  Consider  the  angular  estimate  being  the  search  angle  that 
produces  the  smallest  score  (i.e.,  the  minimum  distance). 
One  way  to  characterize  the  estimation  error  is  to  use  a 
confusion  matrix  in  much  the  same  way  the  confusion 
matrix  is  used  to  characterize  a  classifier. 

•  Provide  a  list  of  discrete  estimates  being  the  search 
angles  for  the  first  few  smallest  scores,  together  with 
their  respective  figure  of  merit.  This  is  an  extension  of 
the  first  estimation  method. 

•  Generate  a  distribution  over  the  search  interval  as  the 
probability  that  each  search  angle  being  the  true  pose. 
This  provides  a  likelihood  function  that  can  be  used  in 
Bayesian  reasoning. 

For  the  first  method,  the  search  angle  that  produces  the 
smallest  score  is  taken  as  the  pose  estimate.  When  the 
index  of  the  lowest  dip  is  searched  over  360  degrees  in 
azimuth  for  the  vehicle  transport,  most  estimates  lie  along 
the  true  azimuth  angle.  However,  there  are  some  large 
excursions  between  the  two  symmetric  axes.  When  the 
search  is  restricted  within  180  degrees  in  azimuth,  thus 
eliminating  the  symmetry,  the  pose  angular  estimates  are 
shown  in  Fig.  7.  Except  for  infrequent  spikes  of  a  dozen  of 
degrees  due  to  minor  symmetry,  the  rest  stays  within  few 
degrees.  Indeed,  except  for  a  large  spike  of  11  degrees  at 
84  azimuth  angle,  other  errors  are  small  within  2  degrees. 
The  standard  deviation  is  1.0826  degrees  and  0.6918 
degrees  after  the  spike  is  removed. 

For  the  generic  tank,  the  discriminant  output  at  six 
consecutive  azimuth  angles  are  shown  in  Fig.  8  where  the 
true  target  azimuth  angle  varies  from  45  to  50  degrees.  The 
discriminant  values  are  the  lowest  around  the  true  azimuth 
angle.  Fig.  9  is  the  shifted  and  normalized  values  around 
the  lowest  dip.  Again,  not  all  dips  locate  at  zero.  These 
normalized  discrimination  curves  do  not  have  the  same 
slopes  and  the  left  and  right  slopes  are  not  symmetric. 

Again,  when  the  index  of  the  lowest  dip  is  searched  over 
360  degrees  in  azimuth,  there  are  large  excursions  between 
the  two  symmetric  axes.  The  pose  angular  estimates  are 
shown  in  Fig.  10  for  the  search  restricted  within  180 
degrees  in  azimuth.  There  is  a  large  excursion  at  120 
degrees.  The  rest  stays  within  few  degrees  of  the  true 
angles.  Indeed,  there  are  a  large  spike  of  46  degrees  at  the 
azimuth  angle  of  120  degrees  and  another  spike  of  13 
degrees  at  the  azimuth  angle  of  100  degrees.  Other  errors 
are  small  within  few  degrees.  The  standard  deviation  is 
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3.8573  degrees,  1.7656  degrees  after  the  major  spike  is 
removed,  and  1.4193  degrees  if  both  spikes  are  removed. 

The  spikes  are  caused  by  the  regions  near  broadside  where 
minor  symmetries  are  present.  Although  the  correct 
template  has  a  low  determinant,  another  template  just 
happens  to  be  slightly  lower,  thus  creating  a  pose  angular 
estimation  error.  It  is  clear  from  these  figures  that  the  tank 
is  less  discriminatory  than  the  vehicle  transport,  because  it 
has  fewer  features.  The  study  illustrates  the  volatility  of 
using  a  classifier  to  estimate  the  pose  at  some  azimuth 
angles.  As  shown  in  [19],  near  broadside  the  target  may 
not  be  visible  at  all,  making  the  spikes  we  observe  above 
irrelevant. 

3.3  Likelihood  Function  and  Updating 

Denote  the  discriminant  output  by  x  =  [x(i),  i  e  0\T  where 
x(i)  =  MTE(P ,  7})  over  the  search  angular  interval  0. 
Denote  the  likelihood  function  by  y  =  [y(7),  i  e  0\T  where 
y(i)  =  Pr{MTE  =  x\  0=i}.  The  proposed  likelihood  function 
is  given  by: 

y  -  a  exp(-x/  a)  (6a) 

a  =  ^Qxp(-x/ a)  (66) 

e 


where  a  is  a  normalization  factor  and  a  is  a  weight  factor. 

This  function  is  designed  so  as  to  inverse  the  score 
function  (a  peak  rather  a  dip  at  a  likely  azimuth),  stretch  it 
(give  a  larger  value  to  more  likely  azimuth),  and  scale  it 
(normalize  all  into  unity).  The  structure  and  the  factor  a  of 
the  transform  are  design  parameters. 

Consider  the  state  space  being  the  discretized  azimuth 
from  0  to  360  degrees  0k  where  k  is  the  time  index. 
Without  any  specific  information  about  the  evolution  of  0k 
over  time,  the  previous  distribution  can  be  used.  Instead  of 
working  with  the  pose  angle  directly,  we  may  choose  to 
work  with  its  distribution  over  the  state  space.  In  this  case, 
the  time  evolution  is  characterized  by  an  angular  spread 
matrix,  denoted  by  A,  as: 


A  =  [A,] 


A„=Pr{0k+l  =  j\Ok=i} 


(7a) 

(7b) 


I>,=1 

;=°  (7c) 

Each  entry  Ajt  in  the  angular  spread  matrix  A  prescribes  the 
probability  that  the  true  pose  at  time  k+ 1  is  j  given  the  true 
pose  at  time  k  is  i.  The  state  equation  can  be  written  as: 

£jfc+ujfc  =  Azk 


--mOk=0\Xk],  Pr{#£  =  A0\  Xk], 


Pr{0k=A0\Xk}]T 


(8a) 


(8b) 


Xk  ={x5,s<£} 


(8c) 


where  AO  is  the  angular  step  used  to  generate  all  the 
templates  for  search  and  N  is  the  number  of  search  steps 
covering  the  angular  uncertainty. 

Denote  the  confusion  matrix  by  C,  which  defined  as: 


C  =  [C,] 


(9a) 


Cj,  =Pv{6  =  j  \  0  =  i) 


(9b) 


N- 1 

IC,=1 


7=0 


(9c) 


It  prescribes  the  probability  that  the  discriminant  indicates 
an  angular  estimate  j  when  the  true  pose  is  actually  i.  The 
confusion  matrix  is  obtained  by  running  a  Monte  Carlo 
simulation  of  the  MTE  algorithm  over  target  templates.  In 
a  sense,  the  confusion  matrix  captures  the  inherent 
resemblance  or  correlation  of  the  target  HRR  profile  as  a 
function  of  azimuth.  This  is  reflected  in  the  local  shape 
(peaks  and  valleys)  of  the  discriminant.  It  also  provides  a 
performance  figure  of  merit  of  a  classifier  as  its  sensitivity 
to  the  azimuth.  Conventional  classifier  designs  actually 
deemphasize  the  sensitivity  to  azimuth  by  using  an 
averaged  template  over  several  degrees  of  azimuth  for 
instance.  In  contrast,  we  seek  for  a  classifier  that  is 
sensitivity  to  azimuth  in  our  application  of  pose  estimation. 
As  a  result,  the  use  of  confusion  matrix  C  is  intended  to 
pick  up  those  distributions  that  are  spread  out  due  to  non¬ 
orthogonality. 

The  predicted  distribution  of  the  discriminant  is: 

lk\k- 1  =CZk\k-i  nm 


Given  the  likelihood  function  yj,  in  (6a),  applying  the 
Bayes’  formula  and  the  total  probability  theorem,  we 
obtain  a  recursive  algorithm  as: 


diag{yk  }ym  i  _  diag{yk  }CAzk_t 

yT y  yTCAZk  i 

Lk  Lk  ~k~l 


(ii) 


The  derivation  is  given  as  follows.  For  the  true  pose  0k  -  i, 
we  have 


zkd)  =  Pr{0*  =i\Xk}  = 


Mxk,xk_t} 


(12a) 


Pr{&  I  6k  =i,Xk^}Pi{0k  =i\  X^jPrjX^) 
YJPr{xk\0k=j,Xk_1}Pr{0k  =i  IXw}Pr{XH} 

'  (12b) 


YJMxk\ek=j}Pv{0k=i\Xk_1] 

(12c) 


ykQzm- 1(;')  _  yjMCA^m 

Z  (7)z«-i  (7)  Z  yt  (7)t CAtk-i  ](;) 

( 1 2d ) 

Assembling  the  expression  in  (12d)  for  all  elements  leads 
to  the  vector  format  in  (11). 
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Fig.  1  Simulation  Environment  with  STARS  & 
SigPred  (FATSO) 


generic  tank  vs.  vehicle  transport 
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Fig.  2  Discriminator  Output  for  a  Generic  Tank 
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Fig.  3  Discriminator  for  Vehicle  Transport 


Fig.  4  Discriminator  Output  for  Generic  Tank 


Fig.  5  Discriminator  at  Discrete  Azimuth 
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Fig.  6  Discriminator  Output  (Details) 
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Fig.  7  Dip  Index  within  180  Degrees 
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Fig.  8  Discriminator  Output  at  Discrete  Azimuth 


Fig.  9  Discriminator  (Shifted  and  Normalized) 


template  azimuth 

Fig.  10  Dip  Index  within  180  Degrees 


Fig.  1 1  From  Discriminant  to  Likelihood 


Fig.  12  Confusion  Matrix  vs.  Likelihood 


Fig.  11  shows  the  discriminant  outputs  (blue-colored) 
when  the  template  azimuth  is  varied  over  360  degrees  for  a 
true  pose  at  121  degrees.  The  discriminant  is  normalized 
while  the  proposed  likelihood  (red-colored)  with  a  -  8  is 
scaled  with  a  factor  of  3  so  as  to  display  them  in  the  same 
plot.  By  exponential  stretching  and  inversing,  the  dips  of 


the  discriminant  are  transformed  into  peaks  of  the 
likelihood  with  other  values  suppressed  to  around  zero. 

Fig.  12  shows  the  121st  row  of  the  confusion  matrix  (blue- 
colored)  with  the  corresponding  likelihood  (green-colored) 
when  the  template  azimuth  is  varied  from  70  to  180 
degrees.  The  major  peaks  match  quite  well  between  the 
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two  curves  with  some  differences  in  sidelobes.  Noise  and 
clutter  contribute  to  mismatches,  which  may  be  averaged 
out  with  multiple  looks  using  the  recursive  algorithm  given 
in  (11). 

or  a  similar  software  program,  the  azimuth  search  interval 
is  pre-discretized  into  a  fixed  grid.  For  the  fixed  azimuth 
grid  points,  all  templates  are  calculated  off-line.  When  the 
pose  angle  is  considered  as  the  state,  the  state  space  is  also 
discretized  into  a  fixed  grid.  The  Markov  chain,  Bayesian 
rule,  and  total  probability  theorem  are  applied  for  recursive 
estimation. 

The  use  of  a  Markov  chain  model  A  to  represent  the 
dynamics  of  pose  angular  changes  may  be  sufficient  if  we 
have  enough  data  (and  time)  to  carry  out  the  recursive 
computation.  In  this  case  with  high  data  rate,  the  mode 
filter  based  on  discrete-time  point  process  [17]  can  be 
applied.  However,  except  for  a  dedicated  tracking  for  fire 
control,  the  sampling  rate  is  typically  low  for  airborne 
surveillance  applications.  An  alternative  way  to  track 
angular  pose  change  is  to  resort  to  direct  search. 

Given  a  target  profile  in  dB  Pk,  the  MTE  score  is 
calculated  against  a  library  of  N  templates  T  =  [T&  /?  =  0, 
...,  N-l]  as: 

Dk  =[MTE(Pk,Tp\  p  =  0,...,  A-l]r  (13) 

The  discriminant  Dk  is  then  converted  into  the  likelihood 
as  defined  in  (6),  rewritten  as: 

hk  =  LH  (Dk,cc)  (14) 

For  two  consecutive  target  profiles  in  dB  Pk  and  Pk+1,  we 
obtain  their  discriminants  Dk  and  Dk+1  using  (13)  and  their 
likelihoods  Lk  and  Lk+1  using  (14).  It  is  possible  to  combine 
the  two  measurements  using  (11)  directly.  However,  a 
seemly  better  method  is  to  find  the  pose  change  first.  This 
is  done  by  the  following  search: 

CO  =  max{Lrk  Lk+1(a>)} 

aeQ.  (15) 

where  co  is  the  pose  angular  change  in  terms  of  azimuth 
search  steps  and  Q  is  the  set  of  possible  pose  changes.  Lk 
and  Lk+1  are  then  combined  using  (11)  with  A  =  I  and  Zk  = 
Lk  and yjc+i  =Lk+i(a>). 

Similar  filters  based  on  belief  theory  can  also  be  derived 
for  pose  estimation  [2]. 

Consider  a  generic  tank  making  constant  turn  with  a  turn 
rate  of  10  degrees  per  sampling  interval.  Six  samples  are 
available  during  the  turn  of  60  degrees. 

Fig.  13  shows  the  discriminants  (i.e.,  the  MTE  classifier’s 
output)  for  the  azimuth  search  angle  varied  from  0  to  90 
degrees  in  1  degree  increment  when  the  true  target  is  at  an 
azimuth  angle  of  20,  30,  40,  50,  60,  and  70  degrees, 
respectively.  The  dip  for  each  curve  indicates  the  estimated 
target  pose  azimuth  angle. 

Fig.  14  shows  the  transformed  discriminants  or 
likelihoods.  Clearly,  the  likelihoods  in  Fig.  14  provide  a 


3.4  Pose  Angular  Tracking 

Due  to  the  heavy  computation  load  associated  with 
preparing  a  target  template  using  the  SigPred  and  STARS 

better  indication  for  the  pose  angle  than  the  discriminants 
in  Fig.  13. 

To  combine  two  successive  profiles,  the  likelihood  of  each 
new  target  profile  is  circularly  shifted.  The  sum  of 
products  between  the  likelihoods  of  template  profile  and 
the  shifted  target  profile  is  taken  for  each  shift.  The 
optimal  shift  is  the  one  that  produces  the  maximum  for  the 
sums  of  likelihood  products.  This  optimal  shift  when 
scaled  by  the  sampling  interval  is  taken  as  an  estimate  of 
the  underlying  turn  rate.  The  element  by  element 
likelihood  products  of  the  template  profile  and  the 
optimally  shifted  target  profile  are  normalized  to  yield  the 
cumulative  PDF  as  shown  in  Fig.  15.  After  two  and  three 
profile  updates  (combining),  the  cumulative  PDF  tends  to 
become  a  point  mass  with  probability  one. 

Fig.  16  shows  the  estimated  shifts  vs.  the  true  shifts  (a 
negative/positive  value  means  a  left/right  shift  along  the  x- 
axis  or  an  upward/  downward  shift  of  a  column  vector  by 
the  same  amount).  The  average  angular  error  is  0.4 
degrees. 

Two  large  errors  occurred  at  the  3rd  and  5th  shifts,  which 
correspond  to  the  4th  and  6th  target  profiles.  Instead  of 
being  peaked,  these  two  profiles  are  flat  around  the  true 
azimuth  as  shown  in  Fig.  14. 

Instead  of  shifting  a  target  profile,  an  alternative  method  is 
to  shift  the  template  profile.  Furthermore,  the  cumulative 
PDF  may  be  spread  with  the  confusion  matrix.  This  is  to 
prevent  a  build-up  at  one  particular  angular  point. 

Fig.  17  shows  the  cumulative  PDF  obtained  using  the 
template  shifting  (forward  combining)  and  confusion 
matrix  spreading.  It  is  clear  that  after  the  confusion  matrix 
spreading,  the  PDF  does  not  amass  at  a  single  point  as  in 
Fig.  15. 

Fig.  18  shows  the  estimated  shifts  vs.  the  true  shifts.  The 
average  angular  error  is  0.4  degrees.  Looking  at  the 
angular  errors  for  this  case  does  not  show  any  advantage  of 
forward  combing  because  we  have  a  very  good  estimate  at 
the  start. 

Now  consider  a  second  scenario  where  we  start  with  an 
ambiguous  estimate.  Fig.  19  shows  the  discriminants  for 
the  azimuth  search  angles  over  0  to  90  degrees  when  the 
true  target  is  at  50,  60,  70,  80,  90,  and  100  degrees, 
respectively.  The  dip  for  each  curve  indicates  the  estimated 
target  pose  azimuth  angle. 

Fig.  20  shows  the  transformed  discriminants  or 
likelihoods.  Clearly,  the  likelihoods  in  Fig.  20  provide  a 
better  indication  for  the  pose  angle  than  the  discriminants 
in  Fig.  19.  As  in  Fig.  19,  the  1st,  3rd,  and  5th  profile 
discriminant  likelihoods  are  flat,  rather  than  peaked  as  for 
the  2nd  and  4th. 
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Fig.  13  Discriminants 


Fig.  14  Likelihoods 


Fig.  15  Cumulative  PDF  (Backward  Combining) 
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Fig.  16  Angular  Shift  per  Sample  (Turn  Rate) 


Fig.  17  Cumulative  PDF  (Forward  Combining)  Fig.  18  Angular  Shift  per  Sample  (Turn  Rate) 
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Fig.  19  Discriminants 


Fig.  20  Likelihoods 


Fig.  21  Cumulative  PDF  (Backward  Combining) 
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Fig.  22  Angular  Shift  per  Sample  (Turn  Rate) 


Fig.  23  Cumulative  PDF  (Forward  Combining)  Fig.  24  Angular  Shift  per  Sample  (Turn  Rate) 


Fig.  21  shows  the  cumulative  PDF  for  backward 
combining,  and  similar  to  Fig.  15,  it  does  become  a  point 
mass  with  probability  one  eventually. 

Fig.  22  shows  the  estimated  shifts  vs.  the  true  shifts,  which 
develop  a  bias  from  the  start  as  expected.  The  average 
angular  error  is  5  degrees. 


Fig.  23  shows  the  cumulative  PDF  obtained  using  the 
template  shifting  (forward  combining)  and  confusion 
matrix  spreading.  It  is  clear  that  after  the  confusion  matrix 
spreading,  the  PDF  does  not  amass  at  a  single  point  as  in 
Fig.  22. 
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Fig.  24  shows  the  estimated  shifts  vs.  the  true  shifts.  The 
average  angular  error  is  1.5  degrees,  which  is  a  smaller 
angular  bias  in  the  turn  rate  estimation.  This  shows  the 
advantage  of  using  the  forward  combination  over  the 
backward  one. 

4  Conclusions 

In  this  paper,  the  use  of  a  minimum  squared  error-based 
classifier,  namely,  the  MTE  algorithm,  for  angular  pose 
estimation  was  investigated.  The  local  behavior  of  the 
classifier  output  was  analyzed,  leading  to  two  estimation 
methods.  The  first  method  took  the  search  angle  that 
produces  the  smallest  score  as  the  pose  estimate  whereas 
the  second  method  was  based  on  the  concept  of  likelihood 
function  converted  from  the  MTE  classifier  output.  The 
likelihood  function  allowed  for  updating  in  a  probabilistic 
manner,  thus  enabling  angular  tracking  of  targets  that 
undergo  maneuvers.  Simulation  results  were  presented  to 
illustrate  the  operation  and  performance  of  these 
algorithms  and  methods. 

The  RF  signatures  used  to  create  the  range  profiles  via 
STARS  in  this  paper  were  generated  by  SigPred  as  part  of 
FATSO.  As  a  public  domain  database,  the  target  models  of 
SigPred  are  of  low-fidelity,  which  contains  a  lot  of  flat 
plates  without  much  detail.  This  tends  to  make  the 
response  strong  if  viewed  at  the  right  angle  but  much 
weaker  otherwise.  The  specific  results  presented  in  this 
paper  may  not  be  representative  of  a  true  target  but  the 
analysis  remains  valid.  It  is  therefore  interesting  to  apply 
the  methods  presented  in  this  paper  to  similar  target 
models  created  by  a  high-fidelity  simulator  such  as 
Xpatch. 

As  the  preliminary  results  for  proof  of  concept,  this  paper 
considered  favorable  SNR  conditions  without  clutter.  More 
realistic  scenarios  with  clutter  are  considered  in  [18]  where 
targets  near  the  clutter  ridge  need  special  attention.  In 
addition  to  the  development  of  pose-sensitive  classifiers,  it 
is  of  interest  to  investigate  pose  estimators  based  on  robust 
nonlinear  mappings  from  range  profiles  to  pose  angles. 
Once  properly  trained,  a  nonlinear  mapping  can  be 
computationally  more  efficient  than  a  template  search  and 
matching  algorithm.  Our  ongoing  effort  is  directed  to  these 
areas. 
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